Unified description of inelastic propensity rules 
for electron transport through nanoscale junctions 
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We present a method to analyze the results of first-principles based calculations of electronic cur- 
rents including inelastic electron-phonon effects. This method allows us to determine the electronic 
and vibrational symmeties in play, and hence to obtain the so-called propensity rules for the studied 
systems. We show that only a few scattering states — namely those belonging to the most transmit- 
ting eigenchannels — need to be considered for a complete description of the electron transport. We 
apply the method on first-principles calculations of four different systems and obtain the propensity 
rules in each case. 



Electronic transport through atomic-size junctions is of 
immense scientific and technological interest. The impor- 
tance of inelastic effects in electronic currents have been 
revealed in several ground-breaking experiments leading 
to the detection and identification of single molecules [1], 
chemical reactions [2, 3], the detection of vibrations in 
atomic wires , the detection of inelastic effects by flu- 
orescence the modification of electron transport in 
nanotubes [6[, the molecular motion induced by elec- 
tronic currents , and the hydrogen detection in atomic 
wires [sj, just to cite a few examples. Of particular im- 
portance due to its spreading use is the case of vibra- 
tional spectroscopy where the conductance changes due 
to phonon emission is measured [13, [HI- This is of- 
ten referred to as point contact spectroscopy or inelastic 
electron tunneling spectroscopy (lETS) [H, However, ex- 
periments alone are not able to give direct insight into the 
fundamental question on how the detailed atomic struc- 
ture correlate with the electrical transport properties. 
There is experimental evidence of approximate selection 
rules {propensity rules [l2[) such that only a small num- 
ber out of the many possible vibrational modes give an 
inelastic signal. These propensity rules yield clues to the 
geometric and electronic structure of the junctions. It is 
therefore of fundamental interest to compare the experi- 
mental results with first-principles calculations. 

Existing calculations of inelastic effects in electron 
transport have been developed either for particular 
cases 
els 



[12I, 13 , 14 1 or for simplified (one- level) mod- 



I5I , lit - First-principles methods capable of treat- 



ing both weak and strong coupling to the electrodes has 
also been developed [3, E,^- However, the results of 
such detailed calculations involve many electronic states 
and vibrational modes. An advanced analysis is therefore 
needed in order to provide insight into the propensity 
rules. 



In this paper we propose a method for analysis of the 
inelastic transport based on just a few selected electronic 
scattering states, namely those belonging to the most 
transmitting eigenchannels at the Fermi energy (sf) [g^]. 
These scattering states typically have the largest ampli- 
tude inside the junction and thus account for the major- 
ity of the electron-phonon (e-ph) scattering. To illustrate 
our method of analysis and to develop an understanding 
of the propensity rules we consider four cases: (i) atomic 
gold-wires, and (ii) molecular junctions, as well as scan- 
ning tunneling microscope (STM) setups in the (iii) reso- 
nant, and (iv) non-resonant limits. The propensity rules 
can in these cases be understood from e-ph induced tran- 
sitions between scattering states of a few eigenchannels. 

Inelastic scattering of electrons in a device under bias 
can be modeled using nonequilibrium Green's functions 
(NEGF) [H, [H, In particular, the lowest or- 

der expansion (LOE) of the NEGF equations provides 
a tractable description of phonon scattering in first- 
principles calculations [13] • This approximation assumes 
a weak e-ph coupling (M) and that the electronic struc- 
ture changes slowly over a phonon energy {hco). It is 
therefore not applicable to strong e-ph coupling. In the 
zero-temperature limit, the conductance is 
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which can be divided into the Landauer term, with the 
transmission r (at Sf) times the conductance quantum 
Go, and inelastic corrections in the conductance from 
each vibrational mode. In this formulation we have sep- 
arated the inelastic contribution in a symmetric term, 
with respect to bias, from phonon absorption and emis- 
sion processes, and an asymmetric term Q^^^^. Q^^^^ 
is small in the cases studied here since it is (i) strictly 
zero for symmetric junctions and (ii) negligible in the 
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FIG. 1: Phase diagram for a one- level model (inset) illus- 
trating the sign of the conductance change at the onset of 
phonon emission. At a given asymmetry factor a the elas- 
tic transmission r has an upper bound Tmax (black line), and 
the inelastic conductance change undergoes a sign change at 
Tcrossover = Tmax/2 (red dashcd line). The one- level model is 
available online [2111. 



STM configuration when very close to or very far from 
resonance [17, 21]. For these reasons we ignore Q^^^^ 
in the following discussion. Each mode thus gives a step 
6>(|eV| — hujx) in the conductance at the phonon energy. 
The magnitude and sign of this step is given by the scat- 
tering rate per excess bias 
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where G is the retarded Green's function at £/, A = 
i{G — G^) the spectral function, and rL,R the couplings 
to the leads. 

It is instructive to use a one-level model to get an un- 
derstanding of the sign of the conductance step. The 
one-level model couples a single electronic level to two 
electronic leads and a localized vibration |l7|, |2l|. A 



phase diagram is plotted in Fig. [T] for the ratio of the 
coupling to the two leads a = Fr/Fl and the transmis- 
sion r at Sf. For this model the maximal transmission 
is Tmax = 4a/(l + a)^ corresponding to the on-resonance 
case. The crossover from a decrease to an increase in the 
conductance is given by the 1/2 rule HEQ, i.e., at 
^crossover = Tmax/2. Back-scattcring dominates the high- 
transmission system leading to a decrease of the conduc- 
tance while forward-scattering leads to an increase in the 
low-transmitting case j24] . 

The scattering rate j^^^ may be interpreted as a 
competition between an inelastic process, first term of 
Eq. ([2]), that increases the current and an elastic correc- 
tion, second term of Eq. ([2]), that decreases the current 



3, In the case of low transmission, where only 

the inelastic term needs to be considered, Troisi et al. 
[lit Gagliardi et al. [13] have extensively discussed 
the propensity rules. 

Instead of trying to understand the complex issues of 
the competition between elastic and inelastic parts of the 
conductance, we have found that the phonon emission 
rate provides a simple way to obtain the lETS propensity 
rules. In the LOE approximation the power deposited 
into the phonon system is given by ^17^ 
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The first term describes electron-hole damping of the vi- 
brations that drives the actual occupation nx towards the 
Bose-Einstein equilibrium value ub- The second term de- 
scribes the heating of the phonon system in terms of the 
emission rate [l7| 
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where Al,r = GrL,RG''' are the partial spectral func- 
tions from the two leads. To provide physical insights we 
rewrite the trace in terms of a complete set of scatter- 
ing states |^z,r) from the left (right) leads. This gives 
Eq. dH) in the form of the physically transparent Fermi's 
golden rule (FGR). It is advantageous to choose the basis 
as eigenchannels [20], i.e., the scattering states belonging 
to the largest transmission. Since the e-ph coupling is es- 
sentially local in space, it is sufficient to evaluate Eq. (|4]) 
using only a few of the most transmitting eigenchannels 
while the reflected scattering states can be ignored. For 
the examples described below, only one to three scat- 
tering states are needed to account for over 90% of the 
phonon scattering. 

To illustrate how the phonon emission rate leads to the 
lETS propensity rules we have performed calculations 
on four experimentally realized systems which we believe 
corresponds to the four 'corners' of the phase diagram, 
see Fig. [T] and O The calculations were performed within 
density functional theory (DFT) [sO] using our extension 
of TranSIESTA as described in Ref. [l^] • Broadening by 
temperature and lock-in modulation Vrms were included. 

The lETS, defined as d^J/d^ V(d//dV), is shown 
in Fig. [2] for: (a) Symmetric low-transmission case, 
an oligo-phenyl-ethylene (OPE) molecule symmetrically 
thiol-bonded to the hollow position on Au(lll) leads. 
The temperature and modulation voltage used in the 
calculation were T = 4.2 K and Kms = 8 meV. As we 
have described previously [l^ , the calculated lETS com- 
pare qualitatively with measurements [9]. (b) Symmet- 
ric high-transmission case, a 7-atom Au chain connected 
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FIG. 2: Calculated (black lines) and experimental (blue lines) 
lETS representing the corners of the phase diagram in Fig. [H 
(a) OPE molecule with Au(lll) leads, (b) Au chain connected 
to Au(lOO) leads, (c) O2 molecule on Ag(llO), and (d) CO 
molecule on Cu(lll). In case (c) the Fermi energy (s/) as 
been shifted manually to match the experiment (dashed red 
line). The experimental data originates from Refs. |l4l. lol. [25I. 
[26I ]. For the STM configurations (c) and (d), the calculated 
lETS is compared with a rescaled d^I /dV^ . 



to Au(lOO) leads, T = 4.2 K, Kms = 1 meV. There is 
quantitative agreement with experiments [J,^. (c) On- 
resonance STM configuration with an O2 molecule on 
Ag(llO) surface displaying a decrease in the conductance 
upon phonon emission [26]. The STM tip is modeled by 
a single Ag atom on a Ag(llO) lead laterally displaced by 
1.6 A corresponding to the experimental situation [26[, 
T = 13 K, Kms = 7 meV. The lETS for O2/Ag(110) is 
shown both for the self-consistent DFT calculation (solid 
black) as well as with Sf shifted manually by —0.6 eV 
with respect to the DFT result (dashed red, low bias 
conductance 750 nA/V). (d) Off-resonance STM config- 
uration with a CO molecule adsorbed on a Cu(lll) sur- 
face [25]. The STM tip is modeled by a single Cu atom 
positioned on a Cu(lll) lead, T = 5 K, Vrms = 2 meV. 

We find that the size of the conductance drop given by 
^LOE ^^Yl approximated with the phonon emission rate 
^FGR OPE and Au-wire cases, the first eigen- 

channel gives the majority of the phonon emission rate, 
see Fig. [3l In this case the propensity rules follows from 
the symmetry of the scattering states which are similar 
to the Bloch states of the corresponding periodic systems 
[20] . To illustrate this understanding of the lETS we have 
published Mathematica notebooks showing the propen- 
sity rules for tight-binding models [2l| . Before comparing 
the theoretical STM-IETS to the experimental, we note 
that the typical tip-sample distances of STM are much 
larger than what is feasible computationally with a local- 
ized basis set. In practice, we work with a significantly 
smaller tip-sample distances than in the experimental sit- 
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FIG. 3: Comparison between the FGR scattering rate, cal- 
culated using only the single-most transmitting eigenchannel, 
with the fuh LOE rate for (a) OPE molecule with Au(lll) 
leads (zero-bias conductance of 1.0 x 10^^/(sV)), and (b) 7- 
atom Au wire (4.8 x 10^^/(sV)). The sign of the FGR rate 
was chosen to agree with the LOE rate. 



nation, while still being in the tunneling regime. 

Results for the O2 on Ag(llO) system is shown in 
Fig.[2jc). In this case our first-principles calculation fails 
to describe the conductance decreases observed experi- 
mentally [30] . The DFT-based energy spectrum does not 
have a molecular resonance at 6/, c.f., the on-resonance 
case of Fig. [TJ However, DFT predicts a tt* resonance 
0.6 eV below Sf \27]. By manually adjusting Sf to this 
resonance we manage to capture some of the qualitative 
features of the lETS. Although this result is suggestive, 
we have no conclusive evidence that the experimental de- 
crease is caused by this mechanism. 

The results for CO/Cu(lll) [30] is shown in Fig.[2l^d) 
and Fig. [H The height of the lETS peaks are quali- 
tatively captured, we predict that 13 % of the trans- 
mitted electrons will emit a frustrated rotation phonon 
compared to 8% found experimentally for CO/Cu(001) 
m, [2^ . The eigenchannels necessary to calculate the 
transition rates 7^^^ are shown in Fig. [D^a). We note 
that the primary eigenchannel is of cr- type, i.e., rotation- 
ally symmetric around the tip-molecule direction, while 
the secondary and tertiary are of 7r-type. Comparing 
the scattering states to the phonon modes. Fig. [H^b), 
the propensity rules naturally follows from the symmetry. 
The frustrated translation (FT^^^^) and frustrated rota- 
tion {F'Rx,y) modes (in the surface x?/-plane) can only 
scatter between a tt- and a a-type eigenchannel, because 
the modes have 7r-character with respect to the trans- 
port direction, i.e., s2 ^ tl and 53 ^ tl. Note that 
the si ^ t2, 3 transitions are less important due to the 
weight of the scattering states (Fig. [4]). Similarly, the 
cr-type FT^ and CO stretch modes only scatter between 
eigenchannels of the same type, e.g., si ^ tl. Table [J 



4 



(a) 



Eigenchannel 1 



tl 





t2 



Eigenchannel 2 




(b) FT(x) 
3 meV 



FR(x) 
34 meV 



FIG. 4: Eigenchannels and phonon modes for CO on Cu(lll). 
(a) Scattering states associated with eigenchannel 1 (left col- 
umn, n = 1.9 X 10~^) and eigenchannel 2 (right column, 
T2 = 0.4 X 10~^). The top (bottom) row of scattering states 
originate from the tip (substrate) side. Eigenchannel 3 (not 
shown, Ts = 0.4 x 10~^) closely resemble those of eigenchan- 
nel 2 rotated by 90 degrees, (b) FR(x) and FT(x) vibrational 
modes (degenerate with FR(y) and FT(y)). 

hcv (meV) FGR (LOE) subs. ^ tip % Mode 

236 0.9 (0.8) s^ ^ t^, i = 1, 2, 3 100 CO stretch 
48 0.3 (0.3) s^^t^,^ = 2,3 95 FT(z) 
35 8.2 (8.0) s3 ^ tl 95 FR(y) 
34 8.3 (8.1) s2 ^ tl 95 FR(x) 
3 5.9 (5.8) s3^ tl 92 FT(y) 
3 6.1 (6.0) s2^ tl 92 FT(x) 

TABLE I: Vibrational modes and inelastic scattering rates 
for CO on Cu(lll). The FGR rate Eq. Q (using all eigen- 
channels) and the LOE rate Eq. ([1]) are given in units of 
10^^{sV)~^ (elastic transmission = 130 x 10^° (sV)~^). The 
dominating transitions between the eigenchannel scattering 
states are indicated along with the fraction in which these 
processes contribute to the total FGR rate. 



confirms these observations and gives a quantitative ac- 
count of the e-ph scattering in the CO/Cu(lll) system. 
It shows that 7^^^ closely approximates the 7^^^ and 
that the main part is given by the first three eigenchan- 
nels. 

We have presented a method to analyze inelastic e- 



ph scattering in terms of eigenchannel scattering states. 
The main advance in the context of first-principles trans- 
port calculations is to bring the description into a natural 
framework where underlying symmetries of the propen- 
sity rules can be understood. Through four different ex- 
amples corresponding to different transport regimes, we 
demonstrated that the eigenchannel analysis addresses 
the propensity rules in a unified way. 
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